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The Colours of Satellite Galaxies in Groups and Clusters 
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ABSTRACT 

Current models of galaxy formation predict satellite galaxies in groups and clusters 
that are redder than observed. We investigate the effect on the colours of satellite 
galaxies produced by the ram pressure stripping of their hot gaseous atmospheres as 
the satellites orbit within their parent halo. We incorporate a model of the stripping 
process based on detailed hydrodynamic simulations within the Durham semi-analytic 
model of galaxy formation. The simulations show that the environment in groups and 
clusters is less aggressive than previously assumed. The main uncertainty in the model 
is the treatment of gas expelled by supcrnovae. With reasonable assumptions for the 
stripping of this material, we find that satellite galaxies are able to retain a significant 
fraction of their hot gas for several Gigayears, thereby replenishing their reservoirs 
of cold, star forming gas and remaining blue for a relatively long period of time. 
A bimodal distribution of galaxy colours, similar to that observed in SDSS data, is 
established and the colours of the satellite galaxies are in good agreement with the 
data. In addition, our model naturally accounts for the observed dependence of satellite 
colours on environment, from small groups to high mass clusters. 

Key words: galaxies: clusters: general; galaxies: evolution - galaxies: fundamental 
parameters: colours - galaxies: luminosity function 



1 INTRODUCTION 

Recent multiwavelength imaging w i th th e Sloan Digi- 
tal Sky survey (SDSS; lYork et all |200Ch has convinc- 
ingly demonstrated that the colour-magnitude distribution 



data (TKauffmann et al. Il993l: lLacev et al.lll993l: ICole et all 



1994 IKauffmann et al 



19991: 



Cole et all l200d: iBenson et 



Someryille fc Primackl Il99fll: 



20021. 120031 : iBaueh et ail 



2005l : ICroton et al.1 12006: Bo wer et al]|2009) . The recent ad 



of galaxies is bimodal dStrateva et al 200 ll : lllogg et all 



l2002l : lBlanton et alll2003l : iBaldrv et alj|2006j l Other galaxy 
properties, such as the star formation rates, disk-to- 
bulge ratios, concentrations, stellar surface mass densi- 
ties and gas mass fractions also show bimodal distri- 



dition of feedback from ac tive galactic nuclei (AGN) in the 
semi-analytic models (e.g.jGranato et al.li2004l : I Bower et al 



a.s 
concentrations, stellar surlace mass 
fractions a l so show bimodal 
butions (IKauffmann et all 120031; iBrinchmann et all J2004 



| 2006| : ICattaneo ei~alll2006l : ICroton et aHl2006J : iMenci et al 



l200q : iKang et all I2006T ) has provided an explanation for 



iBalogh et alll2004lHogg et alll2004lKannappanl2004 ). Elu- 
cidating the origin of these distributions requires under- 
standing whether they arise in situ or whether they are 
environmentally driven (i.e., the "nature" or "nurture" di- 
chotomy) . 

While they are not the only means by which theo- 
rists attempt to tackle the problem of galaxy formation, at 
present semi-analytic models offer the best hope for under- 
standing the origin of the bimodality. The level of sophistica- 
tion of these models has grown rapidly in recent years and so 
too has their ability to match a wide variety of observational 
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three well-known, yet puzzling observational results: (1) the 
steep cut-off at the bright end of the galaxy luminosity func- 
tion; (2) the fact that most massive galaxies today tend to 
be dominated by old, red stellar populations; and (3) the 
absence of classical cooling flows in the centres of massive 
X -ray clusters. In ad dition, some of these models (e.g., that 
of lBower et al] 12009 . hereafter B06, which we adopt as the 
baseline model in this paper) are also able to reproduce the 
evolution of the K-band luminosity and the galaxy st ellar 
mass functions out to high redshift (see iBaughl 120061 and 
references therein for a more in-depth discussion of the suc- 
cesses and limitations of semi-analytic models) . 

In terms of galaxy colours, semi-analytic models with 
AGN feedback are able to produce a clear bimodal separa- 
tion of colours at high luminosities, match the slope of the 
red and blue sequences and explain the absence of massive 
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blue galaxies. However, it has recently become clear that 
these models are unable to match the relative numbers of 
faint s atellite galaxies on the red and b lue sequences. Specif- 
ically, IWeinmann et al.1 (|2006bl 'l and iBaldrv et all (|2006l ) 
have shown that the faint satellite galaxies in groups and 
clusters modelled semi-analytically are on average too red in 
comparison to satellite galaxies of s imilar lumi n osity in the 
SDSS sample. Using DEEP2 data, ICoil et al.l (|2008l ) have 
shown that a similar problem may exist at higher redshifts 
(z ~ 0.7 — 0.9), where the semi-analytic models predict a 
stronger clustering of red galaxies than is observed. These 
results indicate that one or more physical processes (either 
internal or environmentally-driven) that affect the ability of 
the galaxies to form stars are not yet accurately accounted 
for in the models. 

In the present study, we concentrate on the effect of 
environmental processes on the colours of galaxies, particu- 
larly on the role that the ram pressure stripping of the hot 
gas from galaxy haloes plays. At present, the semi-analytic 
models outlined above treat the stripping of the hot gas from 
the haloes of galaxies in a crude fashion, by assuming that 
the gaseous halo is completely and instantaneously stripped 
as the system crosses the virial radius of the (more massive) 
parent system and becomes a satellite. A direct consequence 
of this stripping is that the only fuel that is available for star 
formation once the galaxy becomes a satellite is that which 
resided in a cold disk when the galaxy first fell into the halo. 
As a result, the satellite galaxy experiences a sharp decline 
in its star formation rate and its stellar population becomes 
red over time. This process of cutting off the supply of hot 
gas that would otherwise cool and replenish the reservoir 
of cold, star- forming gas, is sometimes referred to as "stran- 
gulation" or "starvation" (iLarson. Tinslev fc CaldweHlll98(J: 



iBalogh. Navarro fc Morris! 120001 ) . It is clear, however, that 
the maximally-efficient stripping assumed by most current 
semi-analytic models is not realistic, especially in cases 
where the mas s of the satellite is comparable to that of the 
parent system IjWang et al.ll2007f ). 

In a very recent study, iKang fc van den Boschl (|2008l ) 
investigate the effect of modelling the stripping of satellites 
by an exponential decay. They adopt a stripping factor that 
is independent of the properties and orbit of the parent and 
satellite haloes and vary this to achieve the best match to 
the observations. However, although this empirical approach 
is a step forward, one expects the efficiency of the stripping 
to depend on the orbit of the satellite and on the struc- 
tural properties of both the satellite and parent systems. 
In contrast, the model we present here aims to implement 
a physically motivated description for the stripping effect, 
taking into account the relevant properties of the system. 

Recently, iMcCarthv et al.l (|2008T ) carried out a large 
suite of high resolution hydrodynamic simulations of the 
stripping of the hot gaseous haloes of galaxies and presented 
a simple, physically-motivated model that describes the sim- 
ulation results remarkably well. These authors concluded 
that, typically, galaxies are able to retain a significant frac- 
tion of their hot haloes for long periods of time following 
virial crossing. The results of these simulations are in quali- 
tative agreement with a recent Chandra X-ray survey of mas- 
sive satellite galaxies in hot clusters by Sun et al. (2007), who 
found that most ma ssive satellites do in deed have detectable 
hot haloes (see also lJeltema et al.ll2008i ). The implication of 



these results is that the stripping of the hot gaseous haloes is 
much less efficient than previously assumed. Consequently, 
some replenishment of the cold star-forming reservoir (via 
cooling of the hot halo) in the satellite galaxies is expected 
to take place, this prolonging star formation and resulting 
in bluer satellite galaxies. 

In the present study, we incorporate the prescription 
for ra m pressure stripping of hot haloes of IMcCarthv et al.l 
(2008) in the Durham semi-analytic code for galaxy forma- 
tion, GALFORM, and show that this improvement in envi- 
ronmental physics brings the colours of satellite galaxies into 
much better agreement with the observational data. 

The layout of the paper is as follows: In Section [2] we 
present the details of the implementation of the new ram 
pressure model in GALFORM. In Section[3]we present results 
for the fraction of blue galaxies in the new model f £|3.1[l . dis- 
cuss the environmental signatures in the colours of galaxies 
f £|3.2fl and present predictions for the red and blue luminos- 
ity functions at redshift z — 0.1 for galaxies of different type 
(322U- In Section [Hwe summarize and discuss our findings. 



2 A SEMI-ANALYTIC MODEL OF GALAXY 
FORMATION WITH RAM PRESSURE 
STRIPPING OF HOT GASEOUS HALOES 

2.1 The B06 version of GALFORM 

Apart from the minor changes described below, we use the 
version of the GALFORM semi-analytic code described in 
B06. This version of GALFORM makes use of halo merger 
histories extracte d from the Millen nium, Si mulation with 
the t echniques of iHellv et al.l (I2003J) (see also |Harker et al.l 
120061 ). The Millennium Simulation (jSpringel et alJuOOa ) was 
carried out by the Virgo Consortium and it is one of the 
largest simulations of the growth of structure in the ACDM 
cosmology to date, containing approximately 10 billion dark 
matter particles in a cubic volume of (500/i _1 Mpc) 3 with 
a particle mass of 8.6 x 10 8 Mq/L Throughout the paper 
we adopt the same cosmological model assumed in the Mil- 
lennium simulation (and in the B06 model) and quote our 
results in terms of the Hubble variable h = i^o/lOO km 1 
Mpc -1 . 

In terms of baryonic physics, the version of GAL- 
F ORM developed b y B06 has the same basic structure as 
in ICole et al.l (J2000|) , but with the addition of an improved 
treatment of gas cooling and a new scheme for AGN feed- 
back. In the present study, most of the basic parameters of 
the code (including, e.g., the efficiencies of supernovae and 
AGN feedback and the timescales for star formation and 
dynamical friction) are the same as those adopted in B06. 
However, in order to achieve better agreement with the zero- 
point colours of the observed red and blue sequences, we in- 
crease the value of the yield to p = 0.04, which is a factor of 
two higher than the "standard" solar value adopted by B06. 
With this change, the °' (g— r) colours are redder by 0.1 mag 
compared with those in the B06 model (for galaxies on the 
red sequence) . It also improves the metallicity of the intra- 
cluster medium (ICM): whereas in the standard B06 model 
the metallicity of the ICM was too low, Zicm — O.15Z0, in 
the model with double the yields Zicm — 0.3Z© (Bower et 
al. 2008) , a value which is in better agreement with the X-ray 
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measurements (e.g., lBaumgartner et al.ll2005f ). We note that 
large yields such as these have been used in semi-analytical 
models in the past in order to get better agreement between 
the mode l predictions for galaxy colours and the observa- 
tions f see iKauffmann fc Chariot 19981; de Lucia et alj|2004t 
but see also the discussion of lCole et al.ll200d about increas- 
ing yields and consistency with stellar evolution models). 

One further minor change from the version of B06 con- 
cerns the division between haloes in the rapid cooling and 
hydrostatic regimes. B06 adopted a sharp transition between 
these two regimes. In particular, the B06 model assumes 
that radio-mode feedback is only effective in hydrostatic 
haloes with 



£cool(?"cool) > Q cool i ff( 7 "cool), 



(1) 



where £ coo i and r coo i are the cooling time and radius, tg is the 
free- fall time (as defined bv lCole et alJl200Ch and a coo i is an 
adjustable parameter (set in this model to 0.7) that rescales 
the freefall time of the halo and has the effect of controlling 
the position of the break in the luminosity function^. If the 
above condition is satisfied, it is then determined whether or 
not the central AGN is able to inject sufficient power to offset 
the energy being radiated away in the cooling flow and, if the 
available AGN power is greater than the cooling luminosity, 
the cooling flow is assumed to be completely quenched. 

This strict dichotomy in halo properties has the unde- 
sirable effect that small changes in, for example, the gas 
mass fraction of haloes near the transition can greatly affect 
the star formation rate in their central galaxy. This leads to 
a small population of rapidly forming galaxies at the bright 
end of the blue sequence. In the present paper, we have 
refined the criterion to make the transition more gradual. 
To do this, we reduce the cooling rate in haloes for which 
the cooling radius and effective freefall radius [r coo i(i) and 
rs,es{t) = rg (acooli)] are close to each other. Specifically, if 



kcooi - fff.effl < 0.5e cool r coo i, 

the net cooling rate is reduced: 

n-icooLeff = m coo i (0.5 + e coo i[l - r ff /r coo i]) 



(2) 



(3) 



where rh coo i is the cooling rate in the absence of AGN feed- 
back. This function leaves the cooling rate unchanged if the 
inequality above is not satisfied. For large values of e coo i the 
inequality applies for only a narrow range of (r coo i — rfr.eff), 
and the behaviour of the BOB model is maintained. Here we 
adopt e coo i = 10, so that the transition is still quite sharp, 
and the luminosity function is little affected. The effect on 
the colours of galaxies at the bright tip of the blue sequence 
is, however, noticeable: by suppressing gas cooling in these 
objects, the tip of the blue sequence becomes redder, tend- 
ing to curl up towards the red sequence. This provides a 
slightly better match to observational data in this part of 
the colour- magnitude diagram (see SJTT}. 

With these basic parameters, we run a model with 
a complete and instantaneous ram pressure stripping of 
the hot gaseous haloes of satellites (the "default" model) 
and a model where the ram pressure stripping of the hot 
haloes of satellites is calculated using the prescription of 



1 Equation (2) of B06 is incorrect; a coo i should be replaced by 
_— 1 



iMcCarthv et all (|2008f > (henceforth called the "hot ram 
pressure model"). 



2.2 Implementation of Ram Pressure Stripping 

Below we give a brief description of the IMcCarthv et al.l 
1)20081 ) ram pressure stripping model and how it is incorpo- 
rated into GALFORM. The McCarth y et al. model is anal - 
ogous to the original formulation of iGunn fc GottJ (|1972T ) 
for the stripping of a face-on cold disk, except that it ap- 
plies to a spherical distribution of hot gas. Specifically, the 
hot gaseous halo of the satellite will be stripped if the ram 
pressure (P ram ) exceeds the satellite's gravitational restoring 
force per unit area (P gra v): 



2 p G'Mt t, S at(r)pgaB, S at(r) 

am — Pgas,p^ S at ^ -fgrav — CV rp 

r 



P — 

± ram — 



(4) 



where /9 g as. P is the gas density of the parent halo, u sa t is 
the velocity of the satellite with respect to this medium, 
Mtot.sat (r) is the total mass of the satellite within radius r 
and /Ogas.sat (r) is the density of the satellite's hot halo at this 
radius. The coefficient a vp is a geometric constant of order 
unity; McCarthy et al. find that a vp » 2 gives the best fit to 
their hydrodynamic simulations. Note that equation Q has 
not introduced any new free parameters into the GALFORM 
semi-analytical model, as a rp has been tuned to match the 
results of hydrodynamic simulations. 

The satellite-centric radius where P ra m = fgrav is re- 
ferred to as the stripping radius and McCarthy et al. assume 
that within this radius the satellite's gaseous halo remains 
intact while all gas exterior to this radius is stripped on ap- 
proximately a sound crossing time. This physically simple 
model has been shown to match the results of their high 
resolution hydrodynamic simulation^ 2 ] remarkably well for a 
wide range of orbits, mass ratios, and structural properties. 

In our implementation of the hot ram pressure model 
we fix the stripping radius by setting the ram pressure to 
its maximum value, which occurs at the pericentre of the 
satellite's orbit, and the gas is stripped at the instant it 
crosses the virial radius. We adopt this simplification as, 
by default, GALFORM does not track the full orbital evol u- 
tion of the satellite galaxies (but see iBenson et al.l |2003M 3 I. 
In reality, the physical size of the stripping radius is a func- 
tion of time, as the magnitude of the ram pressure varies 
along a non-circular orbit. By setting the ram pressure to 
its maximum value, we overestimate the amount of stripping 
that occurs between the time when the satellite crosses the 



2 This includes both Lagrangian smoo th particle hyd rodynamic 
simulations with the GADGET-2 code iSpringelll2005h and Eule- 
rian adaptive mesh refinement simulations with the FLASH code 
llFrvxell et al.ll2000t) . 

3 The Millennium simulation provides full dynamical informa- 
tion for subhaloes until they are tidally disrupted. In view of 
the limited resolution of the simulation, however, we choose to 
calculate the merging timescale of subhaloes using the standard 

IChandrasekharl J 19431 ) formula rather than following the subhalo 
orbits explicitly, as explained in B06. Since we are interested 
mainly in the broad statistical distribution of galaxy colours, our 
results should not be affected by the loss of information about 
individual orbits. 
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virial radius of the parent halo and when it reaches pericen- 
tre for the first time. However, for typical orbits this is not 
unreasonable, since this timescale is generally a small frac- 
tion of the total amount of time that the satellite spends 
in orbit about the parent halo. In addition, since the full 
orbital evolution of the satellites is not followed, we neglect 
the effects of tidal heating and tidal stripping of the satel- 
lite system. In terms of the removal of hot halo gas, however, 
ram pressure stripping is a more efficient mechanism than 
tidal stripping for the vast majority of satellites (i.e., the 
stripping r adius is typically small er than the satellite's tidal 
radius; see iMcCarthv et alj|200sh . 

The pericentres and the velocities at pericentre (both of 
which are required to compute the maximum ram pressure 
along the orbit) of the satellites are calculated by assuming 
the satellites have the same 2-D joint radial (v r ) and tan- 
gential (vg) velocity distrib ution of infalling substructure as 
measured bv lBensonl (|2005r ) from a large suite of Virgo Con- 
sortium cosmological simulations. iBensonl (|2005l ) finds that 
the following functional form describes the 2-D distribution 
well: 

f(v r ,Ve) = aiveexp[-a 2 (v g - a 9 ) 2 - bi(v$)(v r - b 2 {vg)f] , (5) 

where 01, a 2 , ag, b \(vp) and b 2 (v e) are coefficients and func- 
tions tabulated in IBensonl (|2005r ). We assume that this dis- 
tribution is independent of halo mass and redshift. For each 
satellite we randomly sample this distribution, extracting a 
radial and tangential velocity pair which, in turn, allows us 
to calculate the energy and angular momentum (per unit 
mass) of the orbit. We compute the pericentre radius and 
velocity of the satellite at pericentre by assuming that the 
orbital energy and angular momentum are conserved and 
by treating the satellite as a p oint mass orbiting within 
a Navarro-Frenk- White (NFW) (|Navarro et al.lll996l . Il997h 
potential with the same total mass and concentr ation as 
the parent halo (see, e.g.. lBinnev fc Tremainelll987l ). By as- 
suming that the orbital energy and angular momentum are 
conserved, we are ignoring the dynamical friction force act- 
ing on the satellites. However, in cases where the dynamical 
friction force is strong, the satellite will quickly sink to the 
centre of the parent halo and, in any case, this process is 
important for the most massive satellites which are less af- 
fected by ram pressure stripping. 

By default, the version of GALFORM presented in B06 
assumes that the gas density profiles of the hot haloes of 
all systems (including satell i te sy s tems) follow a /3-model 
jCavaliere fc Fusco-Femianol 1 19761 . Il978h . with fi = 2/3 
(which provides a reasonable match to the X-ray surface 
brightness profiles of massive groups and clusters; Jones & 
Forman 1984), a fixed core radius r c — 0.1R v [ r , and a nor- 
malization po that is set to yield the correct mass of hot gas 
within the virial radius. We adopt the same distribution in 
the present study. We note, however, that we have also ex- 
perimented with NFW density profiles for the hot gas (for 
both the satellite and parent systems), and find that the 
resulting fraction of satellite galaxies that are blue is quite 
similar to the case of the j3- model. This is likely to stem from 
the fact that, typically, the hot gas is only stripped down to 
intermediate radii, where the /3-model and the NFW profile 
are similar. 

In addition to a hot tenuous atmosphere, each halo has a 
reservoir of cold gas that developed from the cooling of the 



hot atmosphere and which can potentially form stars. As 
stars are formed, some of this cold gas is re-heated by feed- 
back from, for example, supernovae winds. Prior to strip- 
ping, the re-heated gas is assumed to follow the same spatial 
and thermodynamic distribution as the hot gas in the halo 
(as in B06). We therefore treat the initial stripping of re- 
heated gas in the same fashion as the hot gas, and transfer 
the same fraction (as computed with the ram pressure strip- 
ping algorithm described above) from the re-heated gas of 
the satellite to the re-heated gas reservoir of the parent halo. 
We record the total mass stripped from the halo as Af str i p 
and the fraction of the hot halo that is stripped as /strip- 
Although this procedure is simple, it is worth noting that 
in reality the spatial and thermodynamic properties of the 
re- heated gas could be quite complicated (e.g., multiphase) 
and that it need not be distributed in the same way as the 
hot gas phase. It would be worth revisiting this issue once 
a better treatment of the re-heated gas is incorporated into 
the semi-analytic model (Benson et al., in prep.). 

The cooling rate of the remaining, unstripped gas is cal- 
culated by cooling only the gas within the stripping radius 
and assuming that stripping does not alter the mean den- 
sity of gas within this radius. We implement this by giving 
the satellite a nominal hot gas mass M^ ot = M^ ot + M atr ip 
(where M^ot is the true hot gas content of the halo) and 
applying the same cooling algorithm as that used for cen- 
tral galaxies (except limiting the maximum cooling radius 
to r-strip rather than i? v ir). This step ensures self-consistency 
in the treatment of the gas cooling between stripped and un- 
stripped galaxies, and therefore that the colours of satellites 
are predicted correctly. 

So far, we have only been considering the stripping of 
gas which was in the hot or re-heated phases at the time the 
satellite halo reached pericentre in the parent halo for the 
first time. However, as long as star formation continues in 
the satellite galaxy, supernova feedback will continue to re- 
heat gas and eject it from the satellite galaxy. The colours 
of satellite galaxies turn out to be sensitive to how much 
of this "secondary" re-heated gas is stripp ed. The numer- 
ical simulations of (|McCarthv et al.l 120081 ) do not provide 
any direct information about the stripping of this re-heated 
gas, since they only treat the stripping of the initial hot 
halo. One can imagine several possible scenarios for the fate 
of the re-heated gas: we adopt a picture in which the the 
re-heated gas is ballistically ejected as relatively cold mate- 
rial that subsequently mixes or evaporates to become part 
of the hot halo. In isolated galaxies, the re-heated gas is 
put back with the same radial distribution as for the orig- 
inal hot gas halo. In satellite galaxies, we can consider two 
extreme cases for the treatment of the secondary re-heated 
gas. Continuing to apply the initial stripping criterion to the 
satellite galaxy as it orbits suggests that the re-heated gas 
should be stripped by the same factor /strip as it is ejected 
from the galaxy. However, applying this at every timestep is 
too extreme: most of the re-heating occurs during the outer 
part of the galaxy's orbit where the ram-pressure force is 
small. It might therefore be more appropriate to consider 
a second case where little of this re-heated gas is stripped 
by ram-pressure effects. Further numerical simulations are 
required to elucidate which of these scenarios is physically 
more realistic and to determine a suitable parameterization 
of the time averaged stripping rate. At this point, the second 
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(minimal stripping) case seems more appropriate since the 
typical orbital timescale exceeds several Gyr, and is com- 
parable to the timescale for the mass-growth of the parent 
halo (see below). 

To allow for these uncertainties, we adopt a partially 
empirical approach, and model the time-dependence of the 
hot gas mass in the satellite halo after its first pericentre 
passage as: 



Afhot — (1 — Estrip /strip) M coo l 

^reheat 



(6) 



Here Mhot is the mass of hot gas available for cooling 
and M coo i is the rate at which it cools. M re heat is the mass 
of gas which has been re-heated by supernovae but not yet 
returned to the hot phase. In the absence of ram-pressure 
stripping, this re-heated gas is assumed to return to the hot 
phase on a timescale T ro i 1C at = idyn/oreheat, as in B06, where 
the value a re heat = 1-26 was chosen to match the observed 
galaxy luminosity function. The effect of ram-pressure strip- 
ping is described by the the term e s tri P /strip, where /strip is 
the stripping factor already calculated for the initial pericen- 
tre using equation (4), and e s tri P is a new parameter (rep- 
resenting the time averaged stripping rate after the initial 
pericentre) which we adjust to fit the observations. The first 
(maximal) stripping case corresponds to e s tri P = 1 if the or- 
bital timescale is much shorter than the re-heating timescale, 
and the second case (minimal stripping) to e s trip = 0. The 
gas which is stripped from each satellite is added to the hot 
gas component of the parent halo. 

We find that if we take e s tri P = 1, then the colour distri- 
bution of satellite galaxies looks very similar to the default 
model in which there is complete stripping of the initial hot 
gas and of the re-heated gas, with all satellites lying on the 
red sequence. For 0.2 < e s tri P < 1, the satellite colour distri- 



bution looks similar to the case e s tr 



1. For estrip < 0.2, 



a blue sequence appears for the satellites, which results in 
better agreement with the observed colour-magnitude distri- 
bution. Values in the range ^ e s tri P < 0.2 result in similar 
colour-magnitude distributions. For this paper, we adopt the 
value estrip =0.1. 

The stripping of satellites is also affected by the growth 
of the halo in which the satellite is orbiting. If we did not 
allow for the effect of halo growth, a small satellite accreted 
at high redshift would not feel the increasing ram-pressure 
effect as the parent halo grows in mass, perhaps becoming 
a galaxy cluster by the present-day. In order to take this ef- 
fect into account, we assign the satellite galaxy new orbital 
parameters and derive a new stripping factor every time the 
halo doubles in mass compared to the initial stripping event. 
If this factor exceeds that applied previously, additional hot 
and re-heated material is removed from the galaxy. This pro- 
cess tends additionally to suppress on-going star formation 
in the satellites of massive haloes. 

An additional physical effect that may be relevant but 
that is not accounted for in our model is that the (ram + 
thermal) pressure force exerted on the hot halo of the satel- 
lite could, in turn, raise the pressure force exerted on the 
cold gas and potential ly stimulate additional star formation 
|Bekki fc Couchl 120031 ) . This would tend to result in bluer 
colours for the satellite galaxies. The magnitude of this ef- 
fect is presently unclear, however, and needs to be quantified 
either observationally or with the aid of hydrodynamic sim- 
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Figure 1. Comparison of the predicted galaxy luminosity func- 
tions of the default model with the hot ram pressure strip- 
ping model at redshift z = 0. The data points represent ob- 
serva tions from the 2dF Galaxy Redshift Survey I INorberg et al,l 
120021) . Magnitudes are bj (Vega) in both the models and the 
observations. Models use volume limited samples and the small- 
est resolved haloes have luminosities ~ 0.03L* (equivalent to 
Mi — 51ogh = —15.7). Thick lines correspond to all galaxies 
and thin lines to satellite galaxies only. 



ulations that accurately include the effects of cooling, star 
formation, and feedback. 



3 RESULTS 

Since it is plausible that including a better treatment of the 
ram pressure stripping may have a fairly significant impact 
on the observable properties of the galaxies, we first check 
if the predicted galaxy luminosity function is significantly 
altered. 

In Fig. [T] we compare the predicted galaxy luminosity 
function for the defauhfj and hot ram pressure stripping 
models at redshift z — 0. The figure shows that our improved 
treatment of the ram pressure does not significantly alter the 
luminosity function. This is primarily because ram pressure 
stripping preferentially affects satellite galaxies, whereas the 
total luminosity function (satellites + centrals) is dominated 
by central galaxies (see H3.3p . In fact, even the luminosity 
function of satellite galaxies alone is only affected by a small 
amount (see thin lines in Fig. [TJ, however, as we show be- 
low, the distribution of satellite colours at fixed luminosity 
is significantly altered. 



4 Our default model luminosity function is essentially the same 
as that of B06, even after the minor changes outlined in £12.21 
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Figure 2. Colour-magnitude diagram (CMD) of z = 0.1 galaxies in the default model (upper panels) and in the hot ram pressure 
stripping model (lower panels). The three panels in each row, from left to right, represent the CMDs of all galaxies, satellites and central 
galaxies, respectively. Magnitudes are SDSS (AB system) g and 01 r at redshift z = 0.1. The contours are spaced linearly in galaxy 
number density, starti ng from 500 per (500 h~ x Mpc) 3 and increasing in levels of 500 per (500 h^ 1 Mpc) 3 . The solid line represents the 
colour cut adopted bv lWeinmann et al] (|2006aj) in order to separate blue and red galaxies in the SDSS. The dashed line represents an 
alternative colour cut (see text). 



3.1 The blue fraction of galaxies 

Fig. [2] shows the predicted colour-magnitude diagram 
(CMD) of present-day galaxies for the default model (up- 
per panels) and for the hot ram pressure stripping model 
(lower panels). The three panels in each row, from left to 
right, represent the CMDs of all galaxies, satellites only, 
and central galaxies only. For a more convenient compar- 
ison with the SDSS data, the model predictions are output 
at redshift z — 0.1 (which corresponds approximately to the 
median redshift of the SDSS sample) and the niters chosen 
for analysis are 01 g and 01 r in the AB system. 

Fig. [2]shows, as expected, that changing the treatment 
of ram pressure stripping primarily affects the colours of 
satellite galaxies. In particular, the distribution of (g — r) 
colours at a fixed magnitude is broader for the hot ram pres- 
sure stripping model and a larger fraction of the satellites 
now have bluer colours. The physical reason for this be- 
haviour is simply that the retention of some of the hot gas 
in the halo of the satellite allows for the replenishment of 
the cold gas reservoir which, in turn, prolongs star forma- 
tion in the satellite after its accretion onto the parent halo. 
Although the satellite galaxies eventually stop forming new 
stars as they consume their gas reservoir, this process is 
much more protracted in the hot ram pressure model. As a 
result, the satellites cover a wider range of colours, filling in 
the region between the two main colour sequences. The red 
sequence is also much less pronounced in the new model. 



To quantify this change in the colours, we classify galax- 
ies as being eit her 'blue' or 'red' by ad opting the colour 
cut proposed bv lWeinmann et al] (|2006al ) for galaxies in the 
SDSS i.e., 

01 (S-^) = 0.7-0.032 ( 01 A/ r - 5 log h + 16.5) , (7) 

This cut, which is represented in the panels of Fig. [2] 
by solid lines, isolates reasonably well the red and the blue 
sequences in both the default and hot ram pressure models 
(see the top left and bottom left panels). This is similar to 
what is observed in the colour-magnitude diagr am of SDSS 
gala xies (see Figure 1 oflWeinmann et alj|2006al ; also Figure 
1 of IWeinmann et alj|2006bh . For the SDSS data, the cut 
adopted by Weinmann et al. tends to isolate the red sequence 
at its base defined at the high luminosity end (as opposed to 
following the minimum between the red and blue sequences). 
This effect is well reproduced by the solid line in Fig. [3] 

We note, however, that adopting exactly the same 
colour cut in our models as in the observations is not neces- 
sarily the best choice. A drawback of using the Weinmann 
et al. colour cut to define blue and red fractions is that 
the results are quite sensitive to the precise position of the 
red sequence. For this reason, we also employ a second cut 
(indicated by the dashed line in Fig. |5J, which intersects 
the model's bimodal contours roughly at t heir minimum. In 
the fol lowing, we will use by default the IWeinmann et al] 
J2006al ) cut and, for illustrative purposes, show also results 
with the alternative cut. 
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The left hand and middle panels of Fig. [3] show, for 
the default and hot ram pressure models respectively, the 
fraction of blue galaxies per ' M r magnitude for satellites 
(filled red circles) , centrals (filled blue squares) and a ll type s 
(filled green triangles) using the IWeinmann et al.l (|2006al ) 
cut. Also shown is the fraction of blue satellites (open red 
circles) using the alternative colour cut. 

Fig. [3] shows that the fraction of satellites that are blue 
in the hot ram pressure model is about twice as high as in 
the default model. As discussed above, the actual fraction 
of blue galaxies depends on the colour cut adopted; however 
the relative difference in the blue fractions between the two 
models is largely indepen dent of the cut. In partic ular, we 
find that using either the IWeinmann et al.l (|2006al ) or the 
alternative cut (or any in between these two) , the fraction of 
satellite galaxies that are blue in the hot ram pressure model 
is approximately 2 — 2.5 times larger than in the default 
model. 

How do the models fare in comparison to the obser- 
vational data? In the right hand panel, we show t he blue 
fraction of galaxies measured by IWeinmann et al.l (|2006bJ) 
using an SDSS group catalogue (kindly provided by S. Wein- 
mann) . The fraction of satellite galaxies that are blue in the 
observations increases up to approximately 60% at the faint 
magnitude end. Encouragingly, the newly implemented ram 
pressure model yields a fraction of blue satellites that is in 
strikingly good agreement with the SDSS results. The satel- 
lite blue fraction rises also to about 60% at the faint end and 
the overall trend in the blue fraction with luminosity is also 
close to that observed. In contrast, the default model yields 
a fraction of blue satellite galaxies at most « 30% (at faint 
magnitudes) . This is simi lar to the frac t ion of blue satellites 
(~ 20%) r eported for the|Croton et al.l (|2006P ) semi-analytic 
model by IWeinmann et al.l ( 2006a ) (see the bottom right 
hand panel of their figure 4). 

We also note that the fraction of bright central galaxies 
that are blue is in broad agreement with the observations. 
This is largely a consequence of including a prescription for 
AGN feedback into GALFORM which tends to halt star for- 
mation in the most massive galaxies. Finally, we note that 
the semi-analytic model predicts many more faint central 
galaxies than observed in the SDSS group catalogue (inde- 
pendently of how the stripping of the hot gaseous haloes is 
treated). However, the lack of faint central galaxies is a selec- 
tion effect of the group sample in the SDSS: the source cata- 
logue excludes groups for which no membe r is brighter than 
° A M r = -19.5 + 5 log ft, where h = 0.73 (jWeinmann et al.1 
l2006al ). Although we do not attempt to reproduce the mock 
groups catalogue in detail here, including this selection cri- 
terion tends to makes the satellite blue fraction shown in 
Fig. [3] greater by 0.1, further improving the match to the 
observations. 



3.2 Environmental dependence of colours 

Ram pressure stripping is an environmentally-driven pro- 
cess, hence one should expect variations in the fraction of 
blue satellites as a function of parent halo mass. The left 
hand panel of Fig. [4] shows the fraction of blue satellites in 
the hot ram pressure model per ' M r bin, calculated using 
the Weinmann et al. colour cut (eq. [7]), for different parent 
haloes mass ranges. The lower masses, 10 12 < Aftot, parent < 



10 14 Mq, correspond to group- type environments, whereas 
the larger masses, M to t, parent > 10 14 Mq correspond to mas- 
sive clusters. For comparison, in the right hand panel we plot 
th e blue fractions o f satelli tes in the SDSS group catalogue 
of IWeinmann et al.l (|2006bl ) split into the same parent halo 
mass bins. 

Reasonably good agreement with the data is obtained. 
The new model reproduces the overall trend that satellites 
of a given magnitude are bluer if they reside in less massive 
systems. The physical explanation for this behaviour is due 
to the slighthy higher density of the hot gas in more mas- 
sive systems (because lower mass systems, on average, have 
converted a larger fraction of their baryons into stars) and 
to the higher orbital velocities of satellites in more massive 
host haloes. The latter effect is the dominant one, especially 
since the ram pressure scales with the square of the satellite's 
velocity. In particular, if we ignore the weak dependence of 
hot gas density on halo mass, the ram pressure scales simply 
as P r am oc t> 2 at . Typically, v sa t is of order the virial circular 
velocity of the parent halo, implying the ram pressure will 
scale roughly as Pram oc Af tot . As a result, stripping is more 
extensive in more massive parent systems, reducing the fuel 
supply for star formation in satellites. 

The model shows a slight increase in the fraction of blue 
satellite galaxies towards the bright magnitude end, i.e. the 
"bump" at -19 > 01 M r > -21. As Fig. g] shows, this effect 
is more pronounced for bright galaxies residing in low mass 
groups, 10 12 < Aftot, parent < 10 13 M Q , and likely reflects a 
limitation of the way in which the transition between rapid 
cooling and hydrostatic regimes is currently treated in the 
code (see N2.ip . 

The direct comparison between the halo masses in the 
models and those inferred from the observations should be 
treated with some caution. In the case of a theoretical model, 
the true mass of the parent halo is known precisely, however 
in the case of the observational da ta one must make use of a 
mass proxy (e.g.. Eke et al.l 2004[) . In the particular case of 
the SDSS data, IWeinmann et al.l (|2006bl ) use an empirical 
relationship between the total optical luminosity of a sys- 
tem and its mass. The semi-analytic models, however, show 
considerable scatter in the relationship between these quan- 
tities. Ultimately, one would like to construct a mock survey 
from the theoretical models with similar characteristics to 
the SDSS catalogue, and analyze the mock data the same 
way as the real data. This is beyond the scope of the present 
paper and will be addressed in a future study. 

We also note that similar results have been obtained 
from studies that attempt to distinguish low mass from high 
mass environment s using other mass proxies. For example, 
iHogg et al.l (|2004r i find that bluer galaxies in the local Uni- 
verse typically reside in low galaxy density environments, 
whereas redder galaxies tend to live in high galaxy density 
environments. T h is is fu lly consistent with the results of 
IWeinmann et al.l (|2006bf ) and our own model predictions if 
one makes the reasonable assumption that the projected sur- 
face densit y of galaxies increases with increasing parent halo 
mass (e.g.. lGladders et alJl200Tt ). 

Lastly, the colours are expected to depend not only on 
the mass of the parent halo but also on the intrinsic prop- 
erties of the satellite. To help disentangle these two factors, 
we plot in Fig.[5]colour histograms (in °' 1 (g — r)) for satellite 
galaxies in both the default (dashed lines) and the hot ram 
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Figure 3. The fraction of blue galaxies per ° M r magnitude. The left panel shows the default model, the middle panel shows the hot 
ram pressure stripping m odel. For comparison, wc plot in the right panel the fraction of blue galaxies in the SDSS data as derived by 
IWeinmann et al.l 1 120060) . The blue fractions are shown separately for satellites (red circles), central galaxies (blue squares) and all types 
(green triangles). The empty symbols represent the fraction of blue satellites with the alternative colour cut (see the dashed line in Fig. 
[2J|. For better comparison, the black dotted line in the middle panel reproduces the blue fraction of SDSS satellites from the right panel. 
All magnitudes arc k-corrected to redshift 0.1. 
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Figure 4. The fraction of blue satellites per M r magnitude calculated using the Weinmann et al colour cut (eq.[7jl, i n parent haloes of 
differen t mass. The hot ram pressure stripping model is shown in the left panel. Data from the SDSS group catalogue of lWeinmann et al,l 
J2006bh are shown in the right panel. In both panels magnitudes have been k-corrected out to redshift 0.1. 



pressure (solid lines) models, divided into different 01 M r 
magnitude and parent halo mass bins. In the new ram pres- 
sure model, satellites brighter than ' M r ~ — 20 tend to 
be red, independent of the parent halo mass (i.e., environ- 
ment). This is because internal processes, specifically AGN 
feedback, are the dominant mechanism for quenching the 
star formation in these systems. Meanwhile, as expected, 
low mass satellites in the same model (corresponding to 
magnitudes fainter than 01 M r ~ —20) tend to be bluer in 
low mass parent haloes and redder in massive clusters. The 
shapes of the histograms and the overa ll trends are quali- 
tatively similar to observational results (JBalogh et al]|2004t 



IWeinmann et alj|2006al ) — but note that the observations in- 
clude both satellites and centrals. For example, the highest 
luminosity galaxies are always dominated by the narrow red 
peak, regardless of environment, a weak tail of blue galaxies 
only becoming visible in the lowest density regions. Since 
we plot only the satellite galaxies here, a distinct blue se- 
quence is not seen; rather, the satellite galaxies increasingly 
occupy transition colours at low halo masses where star for- 
mation has been partially suppressed but not completely 
extinguished. The results of our model suggest that, at the 
low halo mass end, environmental processes are as impor- 
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Figure 5. Fraction of satellite galaxies (f[ 0,1 (g — r), Mtot}) per 01 (g — r) colour bin, divided into different 01 M r magnitude ranges and 
residing in parent haloes of different total mass. Results for the hot ram pressure model are plotted with solid lines and for the default 
model with dot-dashe d lines. Magnitudes incre ase from left to right and parent halo mass increase from top towards bottom. Vertical 
lines correspond to the IWeinmann et al.l Il2006al) colour cut (solid lines) and to the alternative colour cut described in i|3. If dotted lines). 



tant as the intrinsic physical processes in determining the 
colour of satellite galaxies. 



3.3 Luminosity functions of red and blue galaxies 

at z = 0.1 

As a starting point for future comparisons with SDSS and 
other data, we plot colour luminosity functions (LF) for the 
hot ram pressure model. Fig. [5] shows the red and blue LFs 
in the 01 M r band at redshift z = 0.1, where colours are 
separated using the Weinmann et al. cut. Panels include 
all galaxies (left), satellites (middle) and central galaxies 
(right). Dashed lines represent the default model. Current 
observational results for the luminosity functions are not 
able to separate contributions of satellite and central galax- 
ies accurately over this range of luminosity due to the diffi- 
culty of robustly identifying fa int central galaxies (However, 
see the recent HOD analysis of lBrown et alj|200ct ). The com- 
parison in Fig. [6] highlights the origin of the variations of the 
blue fraction seen in Fig. [4] 

With the exception of the blue/red LF of satellites, 
the differences between the default and the hot ram pres- 
sure model are minor. Since the centrals dominate in to- 
tal number density, total LFs are similar for both mod- 



els (cf., discussion of Fig. [T}. When all galaxies are com- 
bined, the blue galaxy LF has a steeper faint end slope than 
the r ed counterpart, in agreement with observational data 
(e.g.. lBaldrv et al.ll200q ) l 5 l. By splitting the luminosity func- 
tion into central and satellite galaxies, we see that this is 
driven by the rapidly increasing preponderance of blue cen- 
tral galaxies at faint magnitudes. In contrast, the luminosity 
function of central red galaxies drops by about three orders 
of magnitude between °M r of —20 and —17. This effect is 
driven by the AGN feedback in the model, and is indepen- 
dent of the stripping model. 

The properties of satellites in the new model are much 
more dependent on the environmental physics included in 
the model. In the new model, the blue and red LFs reach 
similar values at the faint end (consistent with the results 
in Figure [3] showing that at the faint end there are roughly 



5 However, a study that appeared after our paper was submitted 
indicates that the faint end slopes of the red and blue luminosity 
functions in t he SDSS DR6 data are shallo wer than predicted 
by our model I IMontero-Dorta fc Pradal 120081 ). If confirmed, this 
may suggest that other physical pro cesses unaccounted for in our 
model, e.g. tidal stripping of stars I IHenriques et al.ll2008fl , may 
be responsible for the further flattening of the faint end slopes. 
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Figure 6. Luminosity functions in the 01 M r band at redshift z = 0.1. Panels include all galaxies (left), satellites (middle) and central 
galaxies (right). Solid lines correspond to the hot ram pressure model and dashed lines to the default model. Black lines denote all 
galaxies in the (sub)sample, while red and blue lines denote galaxies that are above and below the 01 (g — r) colour cut in equation [7] 
respectively. The green dotted lines in the middle and right hand panels reproduce the total luminosity function of all galaxies for the 
hot ram pressure model (solid black line in left panel). 



similar numbers of red and blue satellites). The dramatic 
increase in the fraction of faint blue satellites can be seen 
by comparing the solid and dashed blue lines. This is an 
alternative way of presenting the information in Fig. [4] and 
underscores the importance of studying the environmental 
dependence of galaxy properties in order to obtain a com- 
plete picture of galaxy formation and evolution. 



4 SUMMARY & DISCUSSION 

Until now, most current semi-analytic models of galaxy for- 
mation have adopted a crude modelling of the ram pres- 
sure stripping of the hot gaseous halos of satellite galaxies. 
In particular, they typically assume complete and instanta- 
neous stripping of the hot gas halo when the galaxy first falls 
in, without regard for the galaxy's mass, orbit, or structural 
properties. This is at odds with results of hydrodynamic sim- 
ulations, and also with recent X-ray observations of galaxies 
in massive clusters. In the present study, we have improved 
the treatment o f strip ping by implementing the model of 
iMcCarthv et al.l (|2008l ) (which has been shown to match 
simulations of ram pressure stripping to high accuracy) into 
the GALFORM semi-analytic model for galaxy formation. 
Although the initial stripping event does not require us to 
add additional parameters to the model, subsequent strip- 
ping of the gas re-heated from the disk requires additional 
parameterization. We parameterize this process by assum- 
ing that most of the ejecta are retained in the galaxy after 
the first pericentre passage. 

We find that the newly implemented treatment of strip- 
ping leads to a significant improvement in the ability of 
the model to match the colours of satellite galaxies in the 
SDSS. The new model is also able to account for the en- 
vironmental dependence of the colours of satellite galaxies. 
Our results suggest that for the majority of satellite galaxies, 
this environmental process can be as important in modify- 



ing the galaxy colours as intrinsic processes, such as AGN 
or supernovae feedback, which operate within the satellites 
themselves. This finding is in broad agreement with previ- 
ous studies that found that internal processes that quench 
star formation do not seem capable of e xplaining the full 
range of colour and morph ology data (e.g.. lWeinmann et al.l 
12006 J : iBaidrv et aljExxSh . 

Although we have only focused on the stripping of hot 
gaseous halos in the present study, other environmental pro- 
cesses may be relevant as well. We now briefly review some 
of these and conclude that none of them appear to be as im- 
portant as the ram pressure stripping of the hot haloes and 
the feedback already incorporated into our present model. 

In cases where the ram pressure stripping of the hot 
halo is complete, some stripping of th e cold gaseous disks 
may also occur (|Gunn fc GottJ 1 19721; lAbadi et all Il999l ; 
lOuilis et al. 20001) . Howe ver, lOkamoto fc Nagashimal (|2003f ) 
and lLanzoni et al.l (|2005l ) have explored ram pressure strip- 
ping of disks in semi-analytic models and have concluded 
that the effect on the colours and star formation rates 
of satellite galaxies are minimal. This most likely stems 
from the fact that disk stripping is only expected to be 
relevant for a minority of satellite galaxies whose orbits 
take them into the very c entres of massive systems (e.g., 
iBriiggen fc De LucialkoOSt ). 

Thermal evaporation of the hot gaseous haloes ( and/o r 
disks) could also be relevant (jCowie fc Songailal Il977l ). 
but observational studies of bubbles and cold fronts 
in X-ray groups and clusters have placed strong con- 
straints on the efficiency o f conduction, concluding that 
it is strongly suppressed (JMarkevitch fc Vikhlininl 120071 ; 
iMcNamara fc Nulsenl 120071 ). Turbulent stripping via the 
generation of Kelvin-Helmholtz and Rayleigh- Taylor insta- 
bilities at the interface between the hot gaseous halo of 
the satellite and the parent system is possible, but the 
tim escale for this type of stripping is generally quite long 



tim escale lor this type or s 
(see IMcCarthv et al.H2008l ) . 
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Other possibly relevant processes include viscous 
stripping (unfortunately at present the viscosity of the 
hot gas in groups and clusters is poorly constrained; 
iMcNamara fc Nulsenl2007l ). tidal effects on the gas as the re- 
sult of the in teraction with the gravitational potent ial of the 
parent halo (JBvrd fc Valtonenlll990l: lMerritdll983J 1 or with 
the ot her satellites (i.e., mergers (|Mihosl ll995T ) or harass- 
ment (|Moore et aI.lll996T 0. and shock heating of the satel- 
lite's hot gas as it fall s in at transonic velocities. However, 
iMcCarthv et alj (|2008P ) have argued that ram pressure strip- 
ping of the hot gaseous haloes is always more efficient than 
tidal stripping by the parent halo's potential or shock heat- 
ing in cases where the satellite mass is less than about 10% 
of the parent halo's mass. 

In this paper we have concentrated only on comparisons 
with data at low redshifts (z < 0.1). In a future study, we 
intend to test the model at other redshifts and compare the 
redshift evolution of galaxy colours with data from current 
deep surveys. Another important application of the model 
is the study of the clustering of galaxies as a function of 
different physical quantities, such as colour. In addition to 
the large scale dependence driven by the relative importance 
of AGN feedback and ram pressure stripping as a function 
of halo mass, our model also predicts a radial dependence 
within a halo driven by the variation in the strength of strip- 
ping with the galaxy orbit. The colour dependence of small- 
scale clustering will constrain the model and perhaps suggest 
a way to improve our treatment of the interaction between 
galaxies and their environment. 

Galaxy catalogues for this model are 

available for download from this http URL: 



http://galaxy-catalogue.dur.ac.uk:8080/Millennium/| 
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